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Abstract. - We develop a model for a driven cell- or artificial membrane in an electrolyte. The 
system is kept far from equilibrium by the application of a DC electric field or by concentration 
gradients, which causes ions to flow through specific ion-conducting units (representing pumps, 
channels or natural pores). We consider the case of planar geometry and Debye-Hiickel regime, 
and obtain the membrane equation of motion within Stokes hydrodynamics. At steady state, 
the applied field causes an accumulation of charges close to the membrane, which, similarly 
to the equilibrium case, can be described with renormalized membrane tension and bending 
modulus. However, as opposed to the equilibrium situation, we find new terms in the membrane 
equation of motion, which arise specifically in the out-of-equilibrium case. We show that these 
terms lead in certain conditions to instabilities. 



Introduction. ~ Phospholipid membranes are a major constituent of the cell. They cover 
the surface of all organelles and play an important role in many fundamental cellular processes 
such as intracellular transport [1]. At a mesoscopic scale, the equilibrium properties of fluid 
membranes can be accounted for by the so-called Helfrich Hamiltonian, which describes the 
mechanics of the membrane in terms of parameters, such as membrane tension, bending 
modulus and spontaneous curvature. These mechanical parameters strongly depend on the 
electrostatic charge of the phospholipid molecules. The electrostatic contributions to the 
bending modulus and surface tension for a fluid membrane have been obtained by several 
authors [2-4] using free energy calculations for simple geometries (sphere, cylinders, planes). 

Real membranes are non-equilibrium systems. They are in general active in the sense 
that they are constantly maintained out of equilibrium either by active proteins (such as 
ATP-consuming enzymes) inside the membrane or by an energy flow due to external param- 
eters (such as a lipid flux [5]). An example of model active membrane has been proposed 
in Ref. [6]. In these experiments, a giant unilamellar vesicle is activated by the inclusion of 
bacteriorhodopsin pumps, which transfer protons unidirectionally across the membrane by 
undergoing light- activated conformational changes. In the same work, a hydrodynamic the- 
ory has also been developed to calculate the nonequilibrium fluctuations of the membrane, 
induced by the activity of the pumps. In this paper, we study theoretically an electrically 
neutral membrane containing passive ion channels in an electrolyte solution and driven out 
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Fig. 1 - Left: Sketch of a quasi-planar fluid membrane embedded in an electrolyte in the presence 
of an applied potential difference and concentration gradients (not represented). Right: equivalent 
electrical circuit, with the contribution of the two ions fc = 1, 2 in parallel, and for each ion, two 
conductances (Gk for the membrane and Hk for the electrolyte) and an electromotive force Ek in 
series. When the finite thickness of the bilayer is included, there is also a capacitor C in parallel 
which contains contributions from Debye layers and from the membrane. 



of equilibrium by the application of a DC electric field or by ion concentration gradients [7] 
(Fig.^. As in refs. [6,8,9], we use a generalized hydrodynamic description appropriate for low 
Reynolds number, that we couple to the electrokinetic equations [10] within Debye-Hiickel ap- 
proximation. We calculate the undulation fluctuations by expansion around the planar state 
of the membrane. With minor modifications, our method can also be used to calculate the 
electrical contribution to the fluctuations of a membrane containing pumps. This contribu- 
tion, was not taken into account in ref. [6]. We first consider the case of a membrane of zero 
thickness, and then generalize the model to the case of a bilayer with a finite thickness and 
a finite dielectric constant, lower than that of the solvent. The later model is more realistic 
since it contains both capacitive effects and ions transport [11, 12]. 

Membrane oj zero thickness driven by the application of an electrostatic potential difference. 
- In this section, we compute the electrostatic potential (f>, and then use the electrostatic 
forces as a source term in the hydrodynamic equations to obtain the membrane equation of 
motion. We use a linear theory and assume a steady state. The quasi-planar membrane is 
located in the plane z = (Fig.QJ, it is embedded in an electrolyte and carries passive channels 
for two types of monovalent ions. There is an imposed potential difference V across the system 
of length L. The boundary conditions for the potential are therefore 4>{z = ±L/2) = ±.V/2 
and dz4>{z — > 0^) = due to the assumption of a vanishing membrane dielectric constant. 
We denote Ck the concentrations of the two ions, where the index fc is I for the positive ion 
[zi = 1) and is 2 for the negative ion {z2 = —1). Far from the membrane these concentrations 
are fixed, so that c^(z — ±L/2) = 71=*=, where the superscript -I- (resp. — ) denotes z > 
(resp. z < 0). Everywhere in the electrolyte we linearize the concentrations around these 
bulk values. A point on the membrane is characterized by a 2D vector field and for small 
undulations a height function /i(r_L). The inverse Debye-Hiickel length is on the positive 
side (resp. k~ on the negative side) with = 2e'^n^ / eksT . Within a linear theory for the 
concentrations and the potential, the fluxes of the positive ions of bulk diffusion coefficient 
Di and of the negative ions of diffusion coefficient D2 are proportional to the ion chemical 
potential gradients 



-n^D2 [Vpf ~ VV'^i , 



(1) 
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where are normalized concentrations such that — n^(l + p^), and ■(/; — ec^jkoT is 
the normahzed potentiaL The charge density is = en^{p^ ~ pt)^ obeys Poisson's 
equation. The conservation of both types of ions in an incompressible fluid imposes that 
dck/dt — — V- Jfc, and therefore at steady state the ions fluxes are constant. At the membrane 
surface, the normal bulk flux of each ion species is equal to the flux through the membrane. 
Throughout this paper, we assume that the membrane has a constant conductance per unit 
area Gk for the two types of ions. In general, the conductance depends in a non linear way 
on the electrostatic potential difference across the membrane [11, 12], but here as a first step 
we assume linear conductivity. In the case of symmetric concentrations = n~ , the charge 
density is an odd function of z and the electric field an even function of z, but no such 
symmetry is present in the general case when n+ ^ . 

The electrostatic potential and the charge density on the positive (resp. negative) side 
depends on the surface charge on the positive electrode (on the same side) ct"'" (resp. cr~), 
which is defined as cr* = zLedz(j>{z — ±L/2). We find that 



1 I h i2 

\Di D2 

1 f ii -12 



(2) 



where ii = eJi • and 12 — ~eJ2 ■ are the electric currents carried by ion species 1 and 2. 
Note that the electroneutrality is satisfied, since p^{z)dz — —cr^ and j'^^ {z)dz = — cr^. 
The non- vanishing charges of the electrodes compensate the overall charge of the electrolyte. 
The electric currents (for Lk^ ^1) ^'^e given by 

1 + GkLkeT/Dkne^ ' ^ ^ 

where Vk ^ V — Ek, Ek ~ —kBT\og{n'^ / n^) / ezk is the Nernst potential of ion fc = 1,2 and 
n = 2n^n~/(n+ +n^). This equation is consistent with the usual electrical representation of 
ion channels, shown in fig. ^ with the contribution of the two ions in parallel, and for each ion, 
two conductances Gk and Hk and an electromotive force Ek in series [11,12]. Hk corresponds 
to the bulk conductance per unit area of each ion Hk = Dkue^ /ksTL, which can be itself 
decomposed in two conductances in series for each side of the electrolyte. 

We now consider fluctuations to first order in the membrane undulation h, and in a quasi- 
steady state, assuming that the pulsation to « DkK^. The calculations are conveniently 
carried out using Fourier transforms, defined as /(q, z) = d?J^±e^'^''^-^ f{r±, z), where q is 
the transverse wave vector for membrane fluctuations. Rotational symmetry imposes that all 
the fields in our problem depend only onq = |q|. We find that, the first order correction to the 
fluxes given by Eq. ^ vanishes, which indicates that there is no correction to the ion conduc- 
tivities due to membrane fluctuation at this order. More precisely, the first-order steady-state 
solution of our electrokinetic equations is then (j)^^^q,z) = ±K^a'^h{q)e-xp{—K^z)/eK'^ (see 

Fig.EJ and p^^^q,z) = ^K'^^h{q) exp{—K'^z)a'^ /k^, where nf — {q^ + n^'^Yl'^. Having de- 
termined perturbatively the ion concentrations and the electric field, we study the membrane 
undulations imposed by the hydrodynamics of the surrounding electrolyte. The velocity field 
is obtained from the Stokes equation in the presence of a force density f acting on the fluid [13], 

TyV^v - VP -(- f = , (4) 

with the incompressibility condition V ■ v = 0. Away from the membrane, the boundary 
condition is v(rj^, z —>■ 00) = 0. We assume an absence of permeation through the membrane 
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Fig. 2 - Solutions of the electrokinetic equations for a membrane of zero thickness and symmetric 
ions concentrations n~ = = n. The parameters used in this figure are V = 50mV, L = 1/xm, 
Gi = G2 = 10t7~Vm^> Di = D2 = lO'^^crnVs and n = 16.6mM. Top: Left: Normalized potential at 
zeroth order 2<^'''^(2)/V as function of 2z/L. Right: Normalized potential at first order in membrane 
undulation <f)'^^\q, z) for a wavevector q = 10®m~^, the normalization factor is h{q)a4:TT/e. Bottom: 
Left: Normalized force along z at first order fi^\q,z). Right: Normalized force along q at first order 
—if^^ {q, z) (where i is the imaginary number). This is for the same wavevector and the normalization 
factor for the forces is h{q)a^ K^iTr/e. 



and the no-slip condition v{r±,z = 0) = (9f /i(rx, <) e^. Here, the active forces [8,9] only act 
on the ions that are transported across the membrane, i.e. i{r±,z) = p{r±,z)'E{r±,z). The 
zeroth order contribution to this force f(°) is balanced by the pressure gradient of Eq. Q). 
One can then derive the fluid velocity v at the membrane surface by solving perturbatively 
Eq. (01). For q <^ k^, and up to first order in the membrane undulation, we find that 
Vz{q,z = 0) = vf{q,z = 0) + vj{q,z = 0) with 



vfiq, . = 0) = ^ + -y'-, - ^y^-, . (5) 




This normal velocity could be derived from an effective free energy density 

F,ff{q) = [Eq^ + Tq^ + Kq^]h\q)/2 (6) 

with S = 3(cr+VK+ + cr-VK")/2e, K = 3(cr+V'«^^ + (J-'^ and F = -2{a+'^ / k+'^ + 
a^'^ j n^'^)l t. Thus, if nonthermal noise can be neglected, the fluctuation spectrum is = 
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ksT/ ((So + + Tq'^ + (Kq + K)q'^), where we have assumed that the membrane has, in 
addition to the electrostatic contribution, a bare tension Sq and a bare bending modulus Kq. 

Let us now discuss the physical significance of this result: 

(i) The first important point is that this derivation of the fluctuation spectrum is independent 
of the origin of the ion fluxes. In particular, it applies also to the case where the ion fluxes 
are produced internally by active pumps, in the absence of any external potential difference 
V = and concentration gradients ^ Alternatively, it also applies when the fluxes 
are produced by concentration gradients ^ in the absence of any external potential 
difference 1^ = 0, provided that the two ions are different (we must require Gi ^ G2 or 
Di 7^ D2 so that cr=^ ^ 0). 

(ii) In addition to the renormalizations of the membrane tension and bending modulus, a 
term proportional to q'^ in the effective free energy is present. This term does not exist 
for an equilibrium membrane in an electrolyte solution, as the electrostatic potential only 
contains even powers of q. It is due to the hydrodynamic couplings, which induce a non- 
analytic dependence on q. Furthermore, the negative sign of F suggests the possibility of an 
instability of the membrane [14]. To see this, let us first consider a free membrane. There is 
no tension, and Et = Eg + S = 0. The membrane is unstable because the average squared 
undulation is negative at low wave vector. In the case of a vesicle, there is a finite tension, 
which can be seen as a Lagrange multiplier ensuring the conservation of its area. The vesicle 
is unstable if the total tension is low enough to satisfy the condition Ef < T"^ / A{Kq + K). This 
instability occurs at a wave vector qc = —r/2{KQ + K) suggesting a periodic deformation of 
the membrane at this characteristic wave vector. As expected, the instability threshold for 
the membrane tension increases with the ion fluxes or the applied potential difference. 

(iii) The value of the renormalized membrane tension E can also be obtained from the zcroth 
order electrostatic potential, with a simple mechanical argument. Indeed, the surface tension 
is related to the Maxwell stress by [15]: 



Note that in the r.h.s. of Eq[7| the second term cancels a contribution proportional to L in 
the first term, which originates from the existence of a pressure gradient in the fluid. 

(iv) For order-of-magnitude estimates, we obtain with the parameters reported in Fig |21 
E = 3.2 • 10-1^ Jm-2, F = -IQ-^* Jm"! and K = lO-^'^keT . Although the ion flux is 
consequent and typical of ion channels, the moduli E, F and K are very small due to the 
strong dependance of these moduli on , which is here only 2.3nm. As we show below, 
these low values also reflect the fact that sofar, we have neglected the bilayer character of the 
membrane and its flnite capacitance. 

Bilayer of finite thickness and finite dielectric constant. - In this section, we extend the 
calculation to a bilayer of flnite thickness d and dielectric constant Cm < e. There is then an 
electrical coupling between the membrane and the surrounding electrolyte, with a strength 
measured by the parameter t = em/iKde) [3,16]. For equilibrium membranes, the importance 
of this coupling is discussed in refs. [17,18]. For the sake of simplicity, we only discuss the 
case of symmetric electrolytes: n~ — , Di = D2 = D, and Gi = G2 = G because we 
intend to provide the full solution and more details in a longer paper. We denote by the 
internal potential and (j) the electrolyte potential. When i ^ 0, the boundary conditions at 
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the membrane are modified, they are now: 

ea,c/.(o)(z ^ ±d/2) = emdAfn{z ^ ±d/2), , 
0(0) (z ^ ±d/2) = (z ^ ±d/2) , ^ ^ 

where the first equation is the continuity condition for the normal electric displacement and 
the second equation is the continuity condition of the potential. With these assumptions, we 
find that the current-voltage relation given in Eq.Ostill holds for any value of t. Furthermore, 
the electric field for z > d/2 (resp. z < —d/2) is Ei^\z) = — cr/e — tr/e exp (k (z + (i/2)) 
(resp. Ei°\z) = -a/e - a/e exp {k {—z + d/2))), while the internal field is constant and 
equal to = — (cr + a)/{tnde). We have introduced a the surface charge on the positive 
electrode, cr, which is related to the electrical current i = ii = 12 in a way similar to Eq. [21by 
a = —2i/K^D, and a defined by ct = Jp p'^{z)dz. These two surface charges are related by 
a/cr = [DK^et - 2G) /2{2t + l)G. Note that only the diffusion time of the ions within a Debye 
layer enters in a whereas a also contains the RC characteristic time of the membrane. 

Using Eq. Hand the equations above, the surface tension can now be written as the sum 
of an internal contribution = —emdEf^^ = — (cr + d)'^ /tne and an external contribution 
Tiout — {~^^ ~ 4crCT + nda^) /eK. The negative contribution I]i„ is known as the Lippmann 
tension [19] and is usually larger in absolute value than Sowt- Because of S^, the total 
membrane tension Eq + ^hi + '^out can become negative at some critical value of the internal 
field Em, which leads to the instabilities discussed in ref. [20]. As discussed in the previous 
section, the surface tension E can also be recovered by solving Stokes equation at first order in 
the membrane height with the appropriate electrostatic forces as source terms. Care must be 
used here as the usual definition of the (bulk) electrostatic force pE holds only inside a medium 
with a homogeneous dielectric constant, but not at the boundaries between two dielectric 
media, where the jump in dielectric constant is associated with a localized force. These 
difficulties can be resolved by calculating the electrostatic force directly with the Maxwell 
stress tensor [21]. From the solution of Stokes equation, the velocity is obtained everywhere 
in the domain |z| > d/2. By extrapolating this velocity at z = 0, an effective free energy of the 
same form as in Eq. ((HJ is obtained, with the same tension S as calculated from Eq.[7| We also 
obtain the moduli F = crij (8 + n^d^ + And) /2K^e and K which is a complicated quadratic 
function of cr and a. 

We now discuss the limit of small conductance G ^ which is of experimental relevance. In 
this limit, cr = since there is no current in the medium i = 0. The membrane is a capacitor of 
surface charge a = —emEm = Vent/ {\ + 2t). This means that the equivalent circuit is made of 
three planar capacitors in series, one being the membrane (of capacitance per unit area Cm/d) 
and the other two corresponding to the Debye layers on each side (of capacitance en per unit 
area), making up a total capacitance G — eKt/{2t+ 1). Now Si„ — —a'^/tne, T.out — — (T^/Ke, 
and F = 0. In this limit, F = since the system is at equilibrium and therefore the effective 
free energy must be of the Helfrich form. With the values of the parameters used in Fig. |5] 
except for G which we take to be G=0 and d = 5nm, e,„/e = 1/40, we obtain t = 3.3 • 10^'^ 
and Kd = 7.4, S„ = -8.4 • 10"^ Jm^^ Eo^t = -1-0 • lO"'^ Jm^^ T ^ and K = O.OllfcsT. 
If we use instead G — lOil^^/m^ a typical value for ion channels, the order of magnitude of 
the tension and bending modulus are unchanged and a small value of F = 8.2 • 10~^° Jm~^ 
is found. This indicates that the capacitor model with G = is a good approximation to 
calculate the moduli in this case. The importance of capacitive effects is confirmed by the fact 
that the values of the moduli obtained here are much larger than the ones obtained previously 
in the case of zero thickness. Furthermore, by varying the ionic strength in the case where ion 
transport is present with G ^ 0, we have found that the capacitor model holds at high ionic 
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strength but become invalid at low ionic strength where ion transport has a stronger impact 
on the moduli. This will be discussed in a future work. 

Conclusions. - In conclusion, we have analyzed the steady-state fluctuations of a mem- 
brane driven out-of-equilibrium by an applied DC electric field or by concentration gradients. 
One of our most notable results is the presence of a new term proportional to in the 
fluctuation spectrum that we interpret as arising from hydrodynamics couplings. For a free 
membrane or for a vesicle of zero thickness and negligible dielectric constant, we have found 
a negative value for T which gives rise to a finite wavelength instability. For a bilayer, we 
recover a known zero wavevcctor instability, when the surface tension becomes negative. Since 
r is positive in this case, the other instability at finite wavelength is absent for a bilayer. We 
have also analyzed the role of capacitive effects, and found that they lead to a dominant con- 
tribution in the electrostatic part of the surface tension and bending modulus at high ionic 
strength. The analysis presented here suggests directions for future study. In particular, it 
would be interesting to go beyond the linear approximation for the electrostatic potential and 
the ion concentrations, and discuss with more details the implications for the dynamics of 
active membranes containing for instance gated or mechano-sensitive ion channels. 
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